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I.  INTRODUCTION 


During  the  period  covered  by  this  Technical  Report,  considerable  advances  have 
been  made  in  many  areas  of  research  under  the  general  title  of  “Mechanics  in  Material 
Space”.  These  advances  are  reported  under  three  main  headings. 

Under  the  heading  of  “Conservation  laws  obtained  using  Noether’s  First 
Theorem”,  results  for  the  various  non-homogeneous  systems  considered  are  obtained  via 
Noether’s  first  theorem!  1],  combined  with  Lie’s  group  theory[2].  The  essence  of  this 
approach  is  as  outlined  in  Section  3  of  the  proposal  that  lead  to  Grant  No.  90-0195.  This 
methodology  for  constructing  conservation  laws  is  applicable  only  to  Lagrangian  systems, 
and  we  had  applied  it  to  four  different  Lagrangian  systems  under  the  sub-headings  of  plane 
elastostatics,  bars,  beams  and  plates. 

Under  the  heading  of  “Conservation  laws  obtained  using  Neutral  Action  Method”, 
results  for  the  various  dissipative  systems,  as  well  as  for  one  Lagrangian  system,  are 
obtained  via  the  Neutral  Action  method.  This  new  methodology  for  the  construction  of 
conservation  laws  has  been  established  with  the  support  from  Grant  No.  90-0195,  and  this 
achievement  has  been  reported  in  Section  7  of  the  progress  report  submitted  for  the  period 
April  1,  1990  through  March  31,  1991.  Details  on  this  new  methodology,  which  is 
applicable  to  dissipative  as  well  as  to  Lagrangian  systems,  can  be  found  in  the  published 
brief  note[3]  attached  as  an  Appendix.  Using  the  Neutral  Action  method,  results  obtained 
are  reported  under  five  sub-headings:  linear  viscoelasticity,  non-homogeneous  beams, 
systems  under  initial  sness,  Sezawa  beam  and  fluid  mechanics. 

Under  the  heading  of  “Results  on  General  Theory”,  the  results  presented  dealt  with 
the  comparison  of  the  two  methodologies  for  constructing  conservation  laws,  Noether’s 
first  theorem  and  the  Neutral  Action  method.  The  connection  between  the  Neutral  Action 
method  and  the  symmetries  of  the  governing  equations  for  any  system  of  interest  was  also 
established. 
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Until  the  discovery  of  the  Neutral  Action  method,  there  existed  no  systematic 
procedure  for  constructing  conservation  laws  valid  for  dissipative  systems.  Since 
conservation  laws  are  of  great  value  in  the  analysis  of  fracture  and  defects  in  these  systems, 
most  of  our  effons  during  this  reporting  year  have  been  concentrated  in  developing 
conservation  laws  for  dissipative  systems  using  this  new  methodology.  Results  in  this  area 
are  presented  in  Section  III. 

Results  presented  in  Section  II  and  Section  IV  are  follow-ups  on  the  results 
reported  in  the  progress  report  for  1990-1991.  Any  advances  in  the  areas  covered  within 
these  sections,  as  well  as  the  current  status  of  the  relevant  manuscripts  are  described  in 
these  sec ’  ns. 

During  the  reporting  period  a  Ph.D.  dissertation  was  completed  under  the  title, 
"Conservation  laws  in  non-homogeneous  and  dissipative  systems"  by  Nelly  Y,  Chien, 
who  graduated  in  June  1992. 


II.  CONSERVATION  LAWS  OBTAINED  USING 
NOETHER’S  FIRST  THEOREM 


1)  Conservation  laws  (and  path>independent  integrals)  in  non- 
homogeneous  plane  elastostatics 


As  mentioned  in  the  last  progress  report,  previous  research  by  Eshelby[4], 
Sanders[5],  Rice[6],  Gunther[7],  Knowles  and  Stemberg[8]  relating  to  path-independent 
integrals  (J,  L  and  M)  plays  a  prominent  role  in  the  study  of  fracture  mechanics  of 
homogeneous  elastic  bodies. 

The  aforementioned  studies  are  applicable  only  to  homogeneous  elastic  bodies,  in 
order  to  analyze  fracture  and  defects  in  non-homogeneous  bodies  with  continuously 
varying  elastic  properties,  our  study  has  concentrated  in  developing  conservation  laws  in 
such  non-homogeneous  materials.  Results  obtained  can  be  applied  to  a  novel  advanced 
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material  called  functionally  gradient  material  (FGM)[9].  These  materials  are  expected  to  be 
used  in  light-weight  structures  such  as  aircraft  and  to  solve  problems  in  the  thermal- 
protection  systems  of  aerospace  vehicles. 

Our  mathematical  apparatus,  used  in  constructing  conservation  laws  for  elastic 
bodies  with  continuously  varying  material  properties,  relied  on  Noether’s  first  theorem 
combined  with  Lie’s  group  theory.  Using  a  novel  application  of  this  apparatus,  we 
managed  to  relate  the  condition  of  the  invariance  of  the  Lagrangian  in  Noether’s  theorem  to 
a  condition  relating  the  elastic  moduli.  In  particular,  we  extended  the  classical  J-integral  to 
materials  for  which  Poisson’s  ratio  is  constant,  while  the  Young’s  modulus  varies 
exponentially  or  as  a  power  law.  This  J-integral  is  path-independent  and  characterizes  the 
crack  propagation  in  such  materials. 

However,  we  would  like  to  point  out  that  our  work  has  not  been  exhaustive,  since 
we  restricted  it  to  the  so-called  geometric  symmetries.  It  is  expected  that  by  including 
generalized  symmetries  the  class  of  non-homogeneous  materials,  for  which  we  can  derive 
conservation  laws,  will  be  enlarged. 

A  manuscript  [Ml]*  covering  this  work  has  been  written  up  during  the  reporting 
period  and  will  be  submitted  to  an  appropriate  journal  for  publication. 

2)  Conservation  laws  for  non-homogeneous  bars 

As  mentioned  in  the  last  progress  report,  using  Noether’s  first  theorem  combined 
with  Lie’s  group  theory,  conservation  laws  were  obtained  for  smoothly  non-homogeneous 
bars.  The  non-homogeneity  can  be  due  to  variable  cross-section  or  to  bars  formed  from 
non-homogeneous  materials.  These  conservation  laws  will  prove  useful  in  studying 

Numbers  in  bracket  preceded  by  letter  M  refer  to  the  list  of  manuscripts  at  the  end  of  this  report. 


concentrated  defects,  such  as  jump  discontinuities  and  cracks  in  smoothly  non- 
homogeneous  bars. 

A  manuscript  [M2]  covering  this  work  has  been  written  up  during  this  reporting 
period  and  has  been  accepted  for  publication  in  Acta  Mechamca. 

3)  Conservation  laws  for  non-homogeneous  beams 

As  mentionea  in  the  last  progress  report,  conservation  laws  for  smoothly  non- 
homogeneous  Bemoulli-Euler  beams  as  well  as  Timoshenko  beams,  has  been  established 
based  on  Noether’s  first  theorem  combined  with  Lie’s  group  theory.  The  two  beam 
theories  differ  in  that  the  Timoshenko’s  theory  of  flexural  motions  in  an  elastic  beam  takes 
rotatory  inertia  and  transverse-shear  deformation  into  account  and  contains  two  independent 
variables  instead  of  the  one  transverse  displacement  of  the  Bemoulli-Euler  theory.  Both 
statics  and  dynamics  of  these  beam  theories  are  considered  in  the  derivation  of  conservation 
laws. 

Conservation  laws  are  found  to  exist  only  for  special  classes  of  non-homogeneous 
beams.  These  admissible  non-homogeneities  are  expressible  in  the  form  a  differential 
equation.  Given  any  non-homogeneity  of  a  beam,  one  only  needs  to  introduce  the  material 
dependence  into  the  differential  equation  to  determine  the  conservation  laws  applicable,  if 
any. 

As  previously  reported,  a  total  of  six  conservation  laws  are  found  for  the  static  of  a 
Bemoulli-Euler  beam,  and  a  total  of  eight  for  the  dynamics  of  a  Bemoulli-Euler  beam. 
These  laws  are  obtained  assuming  that  the  beam  is  loaded  only  at  its  end. 

During  this  reporting  period,  conservation  laws  are  obtained  for  the  statics  and 
dynamics  of  a  non-homogeneous  Bemoulli-Euler  beam  under  a  continuously  distributed 
loading.  As  the  non-homogeneities  of  a  beam  admitting  conservation  laws  aie  expressible 
in  the  form  of  a  differential  equation,  the  distributed  loading  on  a  beam  that  admits 
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conservation  laws  is  also  found  to  be  expressible  in  the  form  of  a  differential  equation.  A 
total  of  two  conservation  laws  are  obtained  for  the  static  case,  and  .  total  of  four  for  the 
dynamic  case. 

For  the  non-horn ogeneous  Timoshenko  beam,  it  has  been  noted  in  the  last  progress 
repon  that  four  conservation  laws  have  been  obtained  for  the  static  case,  and  six  for  the 
dynamic  case.  During  this  reporting  period,  two  more  conservation  laws  were  obtained  for 
the  static  case,  and  one  additional  for  the  dynamic  case.  Also,  the  class  of  non¬ 
homogeneities  admitting  conservation  laws  has  been  enlarged. 

However,  we  would  like  to  point  out  that  our  work  completed  to  date  has  not  been 
exhaustive,  since  we  restricted  it  also  to  the  so-called  geometric  symmetries.  It  is  expected 
that  by  including  generalized  symmetries  more  conservation  laws  will  be  obtained,  and  the 
class  of  non-homogeneous  beams  admitting  conservation  laws  will  also  be  enlarged. 

A  manuscript  [M3]  covering  these  results  has  been  accepted  for  publication  in  the 
International  Journal  of  Solids  and  Structures. 

4 )  Conservation  laws  for  non-hontogeneous  plates 

As  mentioned  in  the  last  progress  report,  conservation  Iriws  were  obtained  for 
smoothly  non-homogeneous  Mindlin  plate  using  Noether’s  first  theorem  combined  with 
Lie’s  group  theory.  Mindlin ’s  plate  bending  theory  is  one  that  includes  the  effects  of 
transverse-shear  deformation  and  rotatory  inertia.  As  in  the  case  for  non-homogeneous 
beams,  the  non-homogeneity  of  the  plate  admitting  conservation  laws  is  expressible  in  the 
form  of  a  set  of  differential  equations.  Given  any  non-homogeneity  of  the  plate,  and  using 
this  set  of  differential  equations,  the  applicable  conservation  laws,  if  any,  can  be  readily 
determined. 

It  has  been  noted  in  the  last  progress  report  that  a  total  of  four  conservation  laws  are 
found  for  the  statics  of  a  Mindlin  plate,  and  a  total  of  six  for  the  dynamics.  During  this 


reporting  period,  two  additional  conservation  laws  were  obtained  for  the  static  case  and  an 
additional  tliree  were  obtained  for  the  dynamics.  Also,  the  class  of  non-homogeneous 
plates  admitting  conservation  laws  has  been  enlarged.  These  laws  are  obtained  assuming 
that  the  beam  is  loaded  only  at  its  ends. 

Conservation  laws  were  also  obtained  for  the  statics  and  dynamics  of  a  non- 
homogeneous  Mindlin  plate  under  a  continuous  distributed  loading  during  this  reporting 
period.  Similar  to  case  of  a  Bemoulli-Euler  beam  with  loading,  the  distributed  loading  on  a 
plate  is  also  found  to  be  expressible  in  the  form  of  a  differential  equation.  A  total  of  six 
conservation  laws  were  obtained  for  the  statics  case  and  a  total  of  four  for  dynamics. 

Also  established  during  the  reporting  period  were  some  analogues  of  the 
conservation  laws  relating  to  the  J,  L  and  M  integrals  in  elasticity.  These  conservation  laws 
expressed  in  the  form  of  path-independent  integrals  are  obtained  for  the  statics  of  a  Mindlin 
plate  with  and  without  distributed  loading.  For  dynamics  with  and  without  loading,  the 
corresponding  conser\  ation  laws  are  expressible  as  balance  laws,  where  path-integrals  are 
balanced  by  the  rate  of  change  of  volume  integrals.  These  results  are  applicable  in  the 
analysis  of  fracture  and  defects  of  non-homogeneous  plates  similar  to  the  J,  L  and  M 
integral  in  elasticity. 

However,  we  would  like  to  point  out  that  our  work  completed  to  date  has  not  been 
exhaustive,  since  we  restricted  it  again  to  the  so-called  geometric  symmetries.  It  is  expected 
that  by  including  generalized  symmetries  more  conservation  laws  will  be  obtained,  and  the 
class  of  non-homogeneous  plates  admitting  conservation  laws  will  also  be  enlarged. 

A  manuscript  [M4]  covering  these  results  has  been  written  up  during  this  reporting 
period  and  will  be  submitted  to  an  appropriate  journal  for  publication  in  the  near  future. 
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III.  CONSERVATION  LAWS  OBTAINED  USING  THE 
NEUTRAL  ACTION  METHOD 

I)  Conservation  laws  in  linear  viscoelasticity 

As  mentioned  in  the  last  progress  report,  conservation  laws  were  obtained  for  one- 
and  two-dimensional  linear  viscoelasticity  (Voigt  model)  using  the  Neur  al  Action  method. 

A  conservation  law  for  the  one-dimensional  case  reported  previously,  is  that  whose 
time  current  can  be  any  function  of  the  total  stress  and  the  space  current  is  the  negative 
derivative  of  this  function  times  some  derivative  of  the  displacement.  If  one  chooses  the 
unknown  function  as  proportional  to  the  stress  squared,  one  can  obtain  the  results  that 
shows  that  the  dissipation  of  .he  elastic  stress  is  being  equal  to  the  rate  of  work  done  due  to 
tractions  minus  the  energy  dissipation. 

In  addition  to  the  conservation  law  mentioned  above,  two  additional  conservation 
laws  valid  for  one-dimensional  linear  viscoelasticity  were  obtained  during  this  reporting 
period.  One  of  these  laws  is  quite  trivial  in  that  it  expresses  the  fact  that  any  function  of  the 
constant  stress  is  also  a  conserved  quantity.  The  other  conservation  law  expresses  a  relation 
between  stress,  displacement  and  velocity  that  can  be  obtained  by  integrating  the 
constitutive  equation  of  the  syster-  while  taking  into  account  the  equilibrium  equation. 

For  the  two-dimensional  problem,  a  conservation  law  that  relates  the  dissipation  of 
elastic  energy  similar  to  that  in  the  one-dimensional  case  was  previously  obtained  only  for 
two  special  cases,  namely,  a  special  pure  dilatation  case  and  a  second  case  without  the 
dilatation  terms.  However,  during  this  reporting  period,  such  a  law  for  the  complete 
problem  was  finally  obtained. 

In  addition  to  the  conservation  laws  which  describe  the  dissipation  of  energy,  a 
conservation  law  without  a  time  current  is  also  available  for  two-dimensional  linear 
viscoelasticity  as  reported  previously.  Without  a  time  current,  this  law  provides  a  path- 


-8- 


independent  integral  that  might  be  useful  in  the  analysis  of  cracks  b  .d  other  defects  in  linear 
viscoelastic  bodies. 

However,  we  would  like  to  point  out  again  that  our  work  completed  to  date  has  not 
been  exhaustive.  By  allowing  the  characteristic  of  conservation  law  in  the  Neutral  Action 
method  to  depend  on  more  and  higher  order  derivatives  of  the  dependent  variables  of  the 
system,  more  conservation  laws  should  be  obtainable. 

A  manuscript  [M5]  covering  these  results  has  been  "'ritten  up  during  this  reporting 
period  and  accepted  for  publication  in  ZAMP.  A  marked  proof  of  this  paper  is  attached  in 
the  Appendix. 

2)  Conservation  laws  for  non-homogeneous  beams 

I 

During  this  reporting  period,  conservation  laws  were  obtained  for  the  statics  and 

i 

dynamics  of  non-homogeneous  Bemoulli-Euler  beams,  with  and  without  a  distributed 
loading  using  the  Neutral  Action  method. 

Since  the  basic  building  block  for  construttion  of  conservation  laws  by  the  Neutral 
Action  method  are  the  governing  equations  of  the  system  of  interest,  this  method  is 
applicable  to  dissipative  systems  without  a  Lagrangian  function  as  well  as  to  Lagrangian 

i 

systems  governed  by  the  associated  Euler-Lagrange  equations. 

Even  though  the  problem  of  non-homogeneous  Bemoulli-Euler  beams  has  been 
treated  previously  based  on  Noether’s  first  theorem  combined  with  Lie’s  group  theory 
limited  to  geometric  symmetries,  in  order  to  compare  the  two  methodologies  for 
constructing  conservation  laws,  we  applied  the  Neutral  Action  method  to  the  same  system. 

The  resulting  conservation  If  vvs  by  the  Neutral  Action  method  are  shown  to  not 
only  encompass  all  previous  results  derived  for  the  same  system  using  Noether's  first 
theorem  combined  with  Lie’s  group  theory  utilizing  geometric  symmetries,  but  they  are 
also  more  numerous  and  are  applicable  to  a  wider  range  of  non-homogeneities  and  loading. 
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A  total  of  eleven  conservation  laws  are  obtained  for  the  static  case  with  and  without 
loading,  and  a  total  of  fourteen  for  the  dynamic  case  with  and  without  loading. 

The  fact  thai  conservation  laws  deriv-^d  by  the  Neutral  Action  method  encompass  all 
previous  results  obtained  using  Noether’s  first  theorem  combined  with  Lie’s  group  theory 
limited  to  geometric  symmetries,  is  consistent  with  the  comparison  between  the  two 
methodologies.  It  has  been  noted  in  Section  IX  of  the  last  progress  report  that,  for 
Lagrangian  systems,  the  requirement  for  existence  of  conservation  laws  by  the  NA  method 
and  by  Noether's  first  theorem  are  mathematically  identical.  However,  as  this  comparison 
is  based  on  the  general  form  of  Noether's  first  theorem  admitting  an  extension  by  Bessel- 
Hagen[10]  (divergence  symmetries),  it  is  expected  that  the  conservation  laws  obtained  here 
for  non-homogeneous  Bemoulli-Euler  beams  using  the  NA  method  to  be  more  general  and 
to  encompass  all  previous  results  reported  in  Section  ni-3. 

Also,  for  a  conservation  law  obtained,  the  Neutral  Action  method  of  constructing 
conservation  laws  is  shown  to  be  more  efficient  than  the  classical  method  of  Noether’s  first 
theorem  with  its  extension  by  Bessel-Hagen.  While  the  condition  for  existence  of 
conservation  laws  by  the  classical  methods  requires  the  use  of  three  unknown  functions, 
the  Neutral  Action  method  uses  only  one  unknown  function  to  arrive  at  the  same 
conservation  law. 

However,  we  would  like  to  point  out  that  our  work  completed  to  date  has  not  been 
exhaustive.  By  again  allowing  the  characteristic  of  conservation  law  in  the  Neutral  Action 
metliod  to  depend  on  more  and  higher  order  derivatives  of  the  dependent  variables  of  the 
system,  more  conservation  laws  will  be  obtained,  and  the  class  of  non-homogeneous 
beams  admitting  conservation  laws  will  also  be  enlarged. 

A  manuscript  [M3]  covering  these  results  has  been  written  up  and  will  be  published 
in  the  International  Journal  of  Solids  and  Structures. 
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3)  Cons;‘r  i  laws  for  systems  under  initial  stress 

During  this  reporting  period,  conservation  laws  were  sought  using  the  Neutral 
Action  method  for  linear  ftla«itic  systems  under  initial  stress.  Since  structural  components 
are  frequently  prestressed  or  exhibit  some  distribution  of  residual  stresses,  a  conservation 
law  that  relates  to  the  energy  release  rate  will  be  extremely  useful  in  the  analysis  of  fracture 
and  defects  for  such  com.ponents. 

Up  until  the  present,  no  conservation  law  was  yet  derived.  The  difficulties  lie  in  the 
generality  of  the  system  at  hand  (the  initial  stress  state  being  unspecified)  and  the 
determination  of  the  characteristic  functions  for  conscA'ation  laws  which  must  satisfy 
approximately  twenty  partial  differential  equations  simultaneously.  Nevertheless,  it  is 
expected  that  results  will  be  obtained  in  the  near  future  for  some  special  cases  such  as 
systems  under  hydrostatic  or  uniaxial  stress. 

4)  Conservation  laws  for  Sezawa  Beam 

During  this  reporting  period,  conservation  laws  for  the  Sezawa  Beam  were  derived 
using  the  Neutral  Action  method.  This  beam  theory  is  a  mrxiification  of  the  Bemoulli-Euler 
beam  theory  which  incorporates  the  effects  of  internal  damping.  Results  Obtained  indicated 
that  conservation  taws  for  this  system  exist  only  in  physical  space.  As  no  conservation 
laws  were  found  in  materials  space,  a  balance  law  that  relates  to  the  energy  release  rate  is 
unavailable. 

5)  Conservation  laws  for  fluid  mechanics 

During  this  reporting  period,  conservation  laws  were  derived  using  the  Neutral 
Action  method  for  2-dimensional  Navier-Stokes  equations  for  an  incompressible  fluid  with 
constant  fluid  density,  constant  viscosity  and  in  the  absence  of  body  forces.  Results 
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obtained  will  be  most  useful  in  the  numerical  analysis  of  the  system,  as  well  as  in  providing 
insights  into  the  behavior  of  the  fluid,  both  mathematically  and  physically. 

Result  obtained  sc  far  are  still  limited.  In  the  general  case,  a  total  of  four 
conservation  laws  were  obtained.  These  laws  do  not  contain  energy  terms  and  they  relate 
only  to  the  incompressibility  of  the  system.  Also,  conservation  laws  were  obtained  for 
cases  with  special  pressure  fields,  such  as  linear  fields  and  fields  which  satisfy  the  Laplace 
equation.  The  conservation  laws  in  these  special  cases  show  a  balance  of  displacement 
gradients,  pressure  gradients  and  strain  rates.  The  physical  interpretation  of  these  laws  is 
yet  to  be  investigated.  The  search  for  conservation  laws  in  this  area  is  far  from  complete. 
Efforts  will  be  undertaken  to  obtain  additional  conservation  laws  and  to  interpret  their 

I 

physical  significance,  hopefully  on  the  level  of  energies.  1 

I 

IV.  RESULTS  ON  GENERAL  THEORY 
Conservation  laws  and  symmetries 

As  reported  in  the  last  progress  report,  the  Neutral  /iction  method  for  constructing 
conservation  laws  has  a  relation  to  the  concept  of  symmetries.  For  each  characteristic  of  a 
conservation  law  by  the  Neutral  Action  method,  there  exists  in  the  adjoint  field  a 
corresponding  characteristic  for  the  symmetries  of  the  governing  equations  for  the  system 
of  interest.  However,  the  converse  is  not  true. 

Also  mentioned  in  the  last  progress  report  is  that  the  Neutral  Action  method  of 
constructing  conservation  laws,  when  applied  to  Lagrangian  systems,  is  identical  to  the 
method  by  Noether’s  first  theorem  with  Bessel-Hagen’s  extension.  The  conditions  for 
existence  of  conservation  laws  of  both  methodologies  can  be  transformed  into  identical 
form. 

During  this  reporting  period,  it  has  been  further  shown  that,  within  Lagrangian 
systems,  both  conservation  laws  derived  by  the  Neutral  Action  method  and  by  Noether’s 
first  theorem  with  Bessel-Hagcn’s  extension  have  a  one-to-one  correspondence  to  the 
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variational  symmetries  of  the  Lagrangian  functional.  This  is  an  additional  supporting  tact 
for  the  equivalence  of  both  methods  when  applied  to  Lagrangian  systems. 

A  manuscript  [M6]  covering  these  results  has  been  written  up  and  submitted  to  an 


appropriate  journal  for  publication. 
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In  this  brief  note,  a  new  methodologs  is  advanced  with  the  aim  of  establlshine  conservation  laws  for  dissipative  svstems  de¬ 
scribed  bv  partial  dilTerential  equations.  This  methodologv  extends  Noether's  cv.ebraied  llrst  theorem,  which  is  applicable  onlv 
to  svstems  governed  bv  equations  derivable  variationallv  from  a  Lagrangian  Several  simple,  vet  non-trivial.  examples  serve  to 
illustrate  the  proposed  procedure. 


The  purpose  of  this  brief  note  is  to  adv  ance  and 
illustrate  a  procedure  for  constructing  conservation 
laws  (i.e..  divergence-free  expressions)  for  dissipa¬ 
tive  systems  described  by  partial  differential  equa¬ 
tions.  Such  expressions  are  most  useful  for  a  variety 
of  reasons.  Whereas  for  non-dissipative  systems, 
which  arise  from  a  variational  principle.  Noether's 
first  theorem  ( I  j  is  available  to  establish  conser¬ 
vation  laws  in  a  systematic  fashion  (cf.  a  very  com¬ 
plete  presentation  and  discussion  by  Olver  (2] ).  no 
corresponding  methodology  was  known  and  ex¬ 
ploited  for  dissipative  systems,  since  these  might  not 
be  related  to  a  variational  principle  and  thus  Noeth¬ 
er's  procedure  becomes  inapplicable. 

To  keep  this  note  short,  reference  again  has  to  be 
made  to  ref.  [2]  for  all  details. 

In  a  non-dissipative  system,  let  T  be  a  Lagrangian 
and  £(L)=0  the  associated  Euler-Lagrange  differ¬ 
ential  equation.  Based  on  Noether's  theorem  it  can 
be  shown  that  if  one  finds  a  "characteristic"  0  of  a 
conservation  law  ( having  made  use  of  the  underly¬ 
ing  symmetry ).  then  the  associated  conservation  law 
may  be  expressed  as 

Div  P=P\  =  QE(L)  =  0. 

where  P'  (/=1.  2 . n)  is  an  n-tuple  of  functions 

determined  by  Q.  n  is  the  number  of  independent 
variables  and  the  subscript  indicates  total 
differentiation. 


Now  let  a  dissipative  system  be  given  by  the  par¬ 
tial  differential  equation 

A(«)  =  0. 

where  u  is  the  dependent  variable. 

Even  though  this  equation  may  not  be  derivable 
from  a  variational  principle,  still  we  set 

/A(«)  =  £;. 

but  now  fis  considered  not  to  be  pre-determined  as 
Q  was.  but  has  to  be  found  from  the  above  equation. 
We  note  that  since  , /Ain)  is  required  to  be  diver¬ 
gence-free.  it  has  to  be  a  null-Lagrangian.  and  this  is 
the  condition  which  determines  /i  Thus  we  set 

y  =./'A(«) 

and  require  6y  =  0  identically  in  u  in  order  to  find 

The  generalization  of  the  foregoing  procedure  to 
a  system  of  partial  differential  equations  is 
straifehtforward. 

.As  a  first  illustration  consider  the  diffusion 
equation 

A(//)  =  h,-««„  =0 . 

The  condition  b'J  =0  leads  to 

V  =0  . 

which  is.  in  this  instance,  the  adjoint  equation.  The 
associated  conservation  law  is  given  by 
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P' =  aU\ll-flt,)  .  P’  =  fli. 

.As  a  second  illustration  consider  the  nonlinear 
wave  equation 


A  (  »)  =1/;  -t/lt,  =0  . 

The  condition  6  j  =0  leads  here  to 


f,  -  nl\  =  0 

and  to  a  class  of  consersation  laws  with 

e'=rf/lL-^+.v— ), 

V  «+  I  n  J 


where  n  is  any  real  number  except  -2.  -1.0.  and 
C  is  an  arbitrary  constant. 

Lastly,  consider  linear  viscoelasticity  in  the  form 

A («)  =  >'i/„  +»;»>,,  =0  . 

Here  )’  is  Young's  modulus,  t/  the  viscosity  of  a  Voigt 
element.  The  stress  a  is  related  to  the  displacement 
u  by 


a=Yii,+f]u,i . 

One  possible  conservation  law  may  be  obtained  as 


P'=gia) . 

P'  =  -g'  j  u,d.\=-g  {  Yu,  +  i]U„ ) . 

where  g{a)  is  arbitrary. 

If  we  choose  g(a)  =  a-/2Y.  with 


P‘,  +  Pl  =  0 


and  with 

(T‘'=  >■< . 

we  may  obtain 
d  .  ., 

which  has  an  immediate  physical  interpretation  and 
is  a  useful  result. 

It  is  quite  remarkable  that  the  methodology  ad¬ 
vanced  in  this  note  has  not  been  exploited  earlier. 
Indeed,  the  tools  to  establish  our  procedure  arc  all 
available  (cf.  ref.  (2].  proposition  5.33)"  but  they 
have  been  used  to  construct  a  machinery  whose  pur¬ 
pose  is  different  from  the  one  herein  presented. 

It  is  seen  that  the  procedure  to  construct  conser¬ 
vation  laws  advanced  here,  has  led  to  novel  results. 
Further  aspects  of  the  methodology  presented  here 
will  be  discussed  in  forthcoming  papers. 

This  work  was  performed  with  the  support  of  the 
.AF  Office  of  Scientific  Research  and  the  U.S.  De¬ 
partment  of  Energy.  This  support  is  gratefully  ac¬ 
knowledged.  The  authors  would  like  also  to  thank 
Professor  P.  Olver  for  many  valuable  discussions. 

•'  It  is  relevant  to  note  here  that  we  became  aware  of  proposition 
5.33  after  our  work  was  completed. 
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Conservation  laws  for  linear  viscoelasticity 
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1.  Introduction 

Since  the  introduction  of  the  J,  L,  and  M  integrals  in  fracture  mechan¬ 
ics,  the  importance  of  path-independent  integrals  has  been  widely  recog¬ 
nized.  These  path-independent  integrals  are  conservation  laws  of  a 
divergence-free  form  in  material  space.  Classically,  conservation  laws  are 
obtainable  through  Noether’s  first  theorem  [2],  with  extension  of  Bessel- 
Hagcn  (!].  However,  Noether’s  approach  presupposes  the  existence  of  a 
Lagrangian  function  for  the  mechanical  systems  considered,  and  since 
systems  with  damping  and  various  other  dissipative  mechanisms  do  not 
possess  Lagrangian  functions,  Noether’s  theorem  cannot  be  applied  to 
dissipative  systems.  .. 

In  a  recent  brief  note{-titie^“Qn  Conservation  Laws  for  Dissipative 
Systems”  [4],  a  new  approach  for  constructing  conservation  laws  was 
proposed.  Given  a  system  aoverned  bv  a  set  of  differential  eouations.  the 
proposed  procedure,  termed  the[— Neutral  Action  (NA>-Mfrthed^’)  in  [5], 
allows  one  to  systematically  construct  the  divergence-free  quantities  applica¬ 
ble  to  the  system  considered. 

It  is  the  purpose  of  this  present  contribution  to  derive  conservation  laws 
valid  for  one-  and  two-dimensional  linear  viscoelasticity  (Voigt  model) 
using  the  proposed  method.  The  results  presented  are  not  exhaustive,  but  as 
a  limited  set,  these  laws  should  prove  useful  in  the  analysis  of  defects  and 
fracture  in  linear  viscoelastic  material. 


2.  The' 


NeU'VroJt  ^cV!  e  r\ 

(NA)  method 


Given  any  system  with  m  independent  variables  .r'  (i  =  1, 2, . . . ,  m),  n 
dependent  variables  u*  {k  —  [,2, . . .  ,n),  the  governing  equations  can  be 
represented  by 

(1) 
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Using  the  multi-index  notation  introduced  by  Olver  [3],  uj  in  equation 
(I)  represents  all  possible  pth  order  partial  derivatives  of  u*, 


dx^'  dx  •'=•••  dx^'’ ' 


with  J-U\,h . as  an  unordered  p-tuple  of  integers,  \ 

indicating  which  derivatives  arc  being  taken,  and  # 7  =plindicationlhow 
many  derivatives  are  being  taken.  -n,, 

^ny  system  governed  by  a  set  of  differential  equations  as  in  equation 
C^^gutral  Action  (NA)  Method’fproposed  in  14]  states  thar  ir  ic 
possible  to  construct  conservation  laws  valid  for  the  system  in  the  form  of  rr.a.-‘ 


a.-V'voi 


rAi=DiP'=o, 


(2) 


£>(gg).(fa^t  +r77i  ( «ro“Ton'  (3) 

HereL/‘  =£'_lx\  are  the  unknown  characteristics  of  conservation 

laws,  and  £*  is  the  Euler  operator  defined  as 
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OUj 


(4) 


with 


‘vrvStr-f 


representing  all  possible  pth  order  total  derivatives,  and 

i~^)j  =  Dj  for#y  =  even, 

/c: — — —D-i-  for  #7  =  odd. _ _ 

...'V  V_p.^ P*4-,-71S.£«.eiftjiZ<iaovcs5!; _ iAAlcg^  i>\d;ceCVfc  /jN-., 

_^ce  our  objective  IS  to  construct  some  divergence-tree  expressions  out 
and  since  the  Euler  operator  acting  on  any  total  divergence  always 
gives  a  null  result  by^alculus  of  variation,  it  follows  that  equation  (3)  is  a 
jjcqtnrra^nt  for  existince  of  conservation  laws.  Equation  (3)  also  implies 
thatp^'Ay  is  a  null  Lagpq^an  whose  action  integral. 


A  = 


(5) 


^  *  ft  ^  ^ 

has  vanishing  variation  for  any  dependent  variables  u*,  5 A  =  other 


words,  m  order  to  construct  conservation  laws  for  any  system  (dissinative 
or  I^grangian)  governed  by  a  set  of  differential  equations  A'  =  orb'y  the  NA  Tr^. 

method,  we  try  to  construct  a  product  oflT^whose  action  mtegral  doerngr — ^  ' 

change  variationally.  Hence  the  name  pNeutral  Action  Method!  riven  to 
this  procedure.  - - -  ^ 

rvNe-Vkai 
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In  practice,  given  any  set  of  differential  equati^s,  one  only  needs  to 
solve  equation  (3)  for  the  unknown  characteristics^  and  then  proceed  to 
construct  the  conserved  currents  P'  valid  for  the  syst|m  governed  by  this  set 
of  differential  equations.  . 


3.  One  dimensional  linear  viscoelasticity 


Conservation  laws  derived  for  1-D  linear  viscoelasticity  were  given  as  an 
example  in  the  brief  note  “On  Conservation  Law's  for  Dissipative  Systems” 
[4]  to  illustrate  the  construction  of  conservation  laws  via  the  NA  method. 
For  completeness  of  the  present  paper,  these  results  will  be  presented  again 
in  this  section  with  some  additional  details. 

The  governing  equation  for  1-D  linear  viscoelasticity  (Voigt  model)  is 
given  by 


(6) 


>5 


where  Y  is  Young’s  modulus,  /j^the  viscosity  coefneient  of  a  Voigt 
element,  .r  and  f/the  spatial  coordinate  and  time,  and  uiihe  displacement  for 
the  system.  SuEscriptslin  thia  peped  indicate^ifferentiatlon. 

TTie  stress  a  of  this  system  is  related  to  Ihe  displacement  u  by 

Klre. 


Yu, 


(7) 


If  one  (transform  the  variables  in  this  system  such  that 
(t>-Yu-^TlU„  a 


(8) 


the  governing  equation  of  I-D  linear  viscoelasticity  can  be  written  as 
A  =  0„  =  O.  (9) 

Assuming  the  characteristic  of  conservation  laws  for  the  system  to  be 

(10) 


the  condition  for  existence  of  conservation  laws,  equation  (3),  by  the  NA 
method  will  require  that  ^ 

rc/i)  +  (.1) 

Since  the  only  unknown  in  the  above  equation  is  the  characteristic  / 
which  depends  on  <(>,  0.,,  and  it  follows  that  all  coefficients  of  second  and 


cWortyj.  s.  ^ 
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higher  order  derivatives  of  0  in  this  equation  must  be  set  equal  to  zero 
independently.  The  resulting  set  of  equations  is  as  follows; 


blocV^ 


■‘rts'i  s. 


v" 


o\j  tr 


(12) 


After  some  mathematical  manipulatiops,  the  solution  of  the  above  set  of 
equations  is  found  to  be 


(13) 


Since  /'  and  /-  are  arbitrary  functions  of  it  is  permissible  to  rename 
them  as 


(14) 


with  (  )'  denoting  differentiation  with  respect  to  the  argument.  Defining/' 
and  /-  as  such  will  allow  the  conservation  laws  derived  to  have  a  simpler 
appearance. 

With 


/=-r(^x)<^,  +  V(0J  +  C0<^7% 


and 


A  =  (pxx. 


(15) 

(16) 


it  is  possible  to  construct  divergence-free  expressions,  D^P*  +  D,P‘  =  0,  out 
of  the  product  of  /  and  A.  The  resulting  conservation  law  is 


P’  =  '(0x)0.  +  K<P 


P'  =  gi<t>.). 


(17) 
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Since  the  stress,  o  =  (A,,  of  this  system  is  constant  with  respect  to  x,  the 
conservation  law  pertaining  to  only  h{4>x)  ^  Qt 

=  (18) 

implies  that  any  function  of  the  constant  stress  is  a  constant  with  respect  to 
X,  where  A,(i)  is  an  arbitrary  function  of  time. 

\Vith  <T  =  <j>^  =[cons^  the  conservation  law  corresponding  to  C  /O  can 
be  rewritten  as  *^{o<vi.-rou^-r 

?' =  0  -  =  ^j(0.  (19) 

where  /IjiO  is  an  arbitrary  function  in  t.  In  terms  of  stress  and  displacement, 

?*  =  yu  +  ;/u,  - -xff  = /42(0.  (20) 

This  law  expresses  a  relation  between  stress,  displacement  and  [strain-rate]^ 
which  can  be  obtained  by  integrating  the  constitutive  equation  of  the  system 
while  taking  into  account  the  equilibrium  equation. 

Ini  terms  of  stress,  the  divergence-free  expression  that  pertains  to  only 
gi4>x)  5*  0  can  be  writtten  as 


I*  a,  dx, 


where  g  is  any  function  of  stress.  n 

If  one  chooses  g{a)  =a^l2Y,  and  {denotingl  the  elastic  component  of 
s^ess  as  ff*  =(iu],  and  the  strain  ratel^given  b^  e,  -  this  particular 
'conwrvation  law  can  be. written  as 

/  .2\  ^ 

=  «» -»?«?.  (22) 


whose  physical  interpretation  is  that  the  dissipation  of  elastic  energy 
is  equal  to  the  rate  of  work  done  due  to  tractions  {ae,)  minus  the 
energy  dissipation  (;?£?). 

Additional  conservation  laws  for  1-D  linear  viscoelasticity  can  be  ob¬ 
tained  by  assuming  different  dependence  of  stress  in  the  function  g((r). 


,4.  Two  dimensional  linear  viscoelasticity 

Denoting  the  normal  stresses  in  the  x  and  y  direction  as  <r“  and  a",  the 
shear  stress  as  the  equilibrium  condition  for  any  2-D  continuum 
mechanics  problem  is  given  by 

A'=<rr  +  ff7  =  0, 

A*  =s  ffj'  -h  ff J'  =  0. 


(23) 


c/ 


o-f* 


-r*  i 
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_  (Ka-^A  e>><  -VKe  Vo>A~f  rr><xie-T^ 

For  linear  viscoelasticiiyJ'^e  stress  components  arc  “related  to  the 
normal  strains  in  the  x  and  y  direction  (s".  £^‘)  and  the  shear  strain  (£■''')  by 

a“r=  +  2;i)£"  +  Jis'-’  +  (a+  2/?)£"  + 

+  2/r)£'>  +  ;^“  +  (a-i-2^)£r +a£;“  (24) 

2;i£'^ -f  2^£f\ 

where  <1  and  n  are  the  Lame  constants,  at  and  /?  the  viscosity  coefficients  of 
the  system. 

With  u  and  v  being  the  displacements  of  the  system  in  the  x  and  y 
directions,  the  strain  components  are  given  by 

e“:=u,. 

(25) 

Combining  the  above  equations,  the  two  governing  equations  for  2-D 
linear  viscoelasticity  in  terms  of  displacements  are 

A'  =  (;.  +  2/i)u„  +  fiUy^  +  (x  +  pt)v^ 

+  (a  +  2fi)uJ^,  +  Pu,y,  +  (a  +  P)v^y,, 

A^  =  a+2^l)v„+^iv,^  +  iX+^t)u,y 

+  (a  +  2P)Vy%,  +  Pv„,  4-  (a  +  /?)«,„• 

Given  the  governing  equations  for  the  system,  the  condition  for  exis¬ 
tence  of  conser/ation  laws  by  the  NA  method,  as  given  in  equation  (3),  is 

[£(7' A'  -r/^A-)  C-fV  ■(W)S  O.  (27) 


SincTf A"sl  contain  third  derivatives  in  u  and  v,  in  order  to  evaluate  this 
oinduionfSrVrxistence.  one  needs  to  compute  the  third  total  derivative  of 


i  /'*sjlf  one  assumes  some  characteristics  that  depend  on  x,  y,  t,  u,  v  and  also 
on  derivatives  of  u  and  v,  calculating  he  third  total  derivatives  o^7^^s)wouId 
be  virtually  impossible  in  the  absence  of  advanced  computing  devices  with 
large  memory  capacity.  Due  to  this  difficulty  in  evaluating  equation  (27), 
the  general  solution  for  the  characteristics  is  yet  unknown.  However, 
particular  solutions  may  be  found  heuristically. 

By  restricting  the  dependence  ofy'’sAo  x,  y  and  i  only,  the  condition  for 
existence  of  conservation  law  by  tha  NA  method,  equation  (27),  requires 
that  ^  r 

-{■is 


/L  -t-/i,  =  o. 
/L+/;,  =  o, 

n  +A  =1^ 


(28) 


4^' 
T*  S. 
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wherejCj^is  an  arbitrary  func,‘ion  of  t.  The  corresponding  conservation 
law  is 

P*  =/V“ +/V’' -/Ua  +  2/i)«  +  (a  +  2^K] 

-fyUiv  +pv,]  +  Ml  +  otu,] 

+at»,]  -/;[/4I/  +M,]  _ 


P'  =  0., 


-/.^[p«  +  /8u,]  -/liiX  +  2n)v  +  (a  +  2^)^ 


SDacy|which  might  be  useful  in  numerical  analysis  of  cracks  and  defects  for 
*■*  2'D  fecoelastic  material.  p-  -t-o  -VV\g. 

^  ^^^CConservatio^  law  that  relates'^dissioanon  of  elastic  energy  for  2-D  linear 
.  viscoelasticity  can  be  constructed  if  onelconsideredi  a  special  case  where  the 


where  /’  and/^  are  functions  satisfying  equation  (28h 

If /'./*  are  taken  to  be  constant,  the  above  conservation  law  expresses 
the  condition  of  equilibrium  for  the  system.  Due  to  the  absence  of  a  time 
current,  equation  (29)  provides  path-independent  integrals  in  the  material 


A 

* 


c 


-w 


Lame  constants  (A,  /i)  and  viscosity  coefficients  (a,  are  related  by 

lioL—Xfi,  (30) 

which  implies  that 

P  a  Tz  . 

One  solution  for  existence  of  conservation  laws  in  this  special  case  is 
given  by 

/'  =  A[yu,  +  u„], 

r-  4r  ^  1  (^2) 

r  =  Am  +  u,,]. 

and  the  corresponding  conservation  iawjbeing^ 

, - +  u„)a"  +  (yv,  -f-  v„)ff^] 

P”  =  /l((yu,  -I-  v„)a>'’'  -H  (yu,  +  u„)ff'’1 0  cle.\e<r«. 

P’  =  ((«  +  2^((y«,  +  -1-  (yr,  -K  v,.)'-]  (33) 

+  +  Uy,y  +  im  +  Ox.)^ 

'I  +  2at(y»,  +  «^,X7“x  +  «x.)  +  2^(y»,  +  Wx,)(7«»  +  u,,)]. 

[Denoting  the  elastic  energy  of  the  system  as 


(34) 


s 
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where  is  the  strain  energy  density  given  by 

the  consei”ation|law^for  this  special  case  of  =  Xp  can  be  written  as 
^  [  w'']  =  a'-'s"  +  -[(a  T-  +  er’) 


-  Ise-ye^^  -  4^c•7»■^ 


which  expresses  the  balance  between  dissipation  of  elastic  energy,  rate  of 
work  done  due  to  tractions  and  the  total  energy  dissipated  for|3  special  2-D 
viscoelasticity  problem.  ^ 

The  balance  law  for  the  special  2-D  viscoelasticity  problem  where 
pa  =  X^,  equation  (36),  can  be  verified  to  hold  also  in  the  general  case 
without  any  restraint  on  the  Lame  constants  and  viscosity  coefficients.  The 
conservation  law  that  would  yield  this  result  is  found  to  be 

=  A[u,a"  -1-  v,a^’'  +  u<r  “  -f 

P>' =  -r  ^  va^/ +  ua^”] 

P'  =  -A[2fV*  -h  (a  +  2j5)(u;,u„  +  Ujj'.w)  +  a(i;,u„  -i-  u^v,,) 

■  +  +  0,U„  +  UvJ;„)l.  _^<rVK(S<i 

Within  the  framework  of  the  Neutral  ActiorfMetho^,  equation  (37)  can  be 
obtained  if  onTlmo^edl  the  condition  for  existence  of  conservation^lay.^  as 
given  by  equation  (p^to  ^  V 

I  r-/  X I  4  X  I  \  _  A  1  ^  o \ 


given  by  equation  (p'^to 

l£(/W-Hg‘Ai)^  ^ 

In  which  case,  the  corresponding  characteristics  are 

/•  =  «,  ^•  =  «  N  H- 

P  =  y,  g*  =  i>. 


Vexot^ 


5.  Conclusions 

In  this  present  contribution,  conservation  laws  are  presented  for  one- 
and  two-dimensional  linear  viscoelasticity.  Thescdiyergcnce-free  expressions 
are  obtained  using  the  Neutral  Action  (NA)  {^thodl  applicable  for  con¬ 
structing  conservation  laws  for  dissipative  systems  (4].' — .G^e.-tKt)A 

For  1-D  linear  viscoelasticity,  the  non-irivial  conservation  laws  obtained 
took  a  very  general  form.  The  conserved  current  in  time  can  be  any  function 
of  the  total  stress,  and  the  conserved  current  in  space  is  the  negative 
derivative  of  this  function  times  some  derivative  of  the  displacement.  If  one 
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•ary—TT 

choose^h^nknownVunciion  to  be  proportional  to  the  stress  squared,  the  re¬ 
sulting  conservatiohlaw  shows  balance  [or]the  dissipation  of  elastic  energy 
and  the_^train  rate]  of  the  system.  Obviously,  more  conservation  laws  can  be 
“constructed  by  assuming  diffe.ent  dependences  of  stress  in  the  unknown 
function. 

In  2-D  linear  viscoelasticity,  the  conservation  laws  obtained  are  not  as 
general  as  in  the  1-D  case.  For  the  general  a.nd  a  special  case  of  the  2-D 
problem,  conservation  laws  are  constructed  to  show  the  balance  between  the^ 
dissipation  of  energy  and  rate  of  work  done  parallel  to  the  1-D  [ejeampl^ 
Also,  it  is  found  that  conservation  laws  without  the  conserved  time  current 
_daiexist  which(  provide  path-independent  integrals  that  might  be  useful  in 
analysis  ot  cracks  and  other  defects  in  linear  viscoelastic  material. 

The  conservation  laws  derived  here  are  not  e.xhaustive.  If  one  allows  the 
characteristics  of  conservation  laws  to  be  functions  of  higher*dcrivatives  of 
the  dependent  variables,  more  conservation  laws  can  be  obtained.  In  lack  of 
an  advanced  computing  device  with  large  memory  capacity  which  is  neces¬ 
sary  to  solve  for  high  dependence  characteristics,  the  conservation  laws 
presented  here  are  limited.  However,  as  a  limited  set,  these  laws  are 
non-trivial  results  for  viscoelasticity  and  they  serve  as  e.xaraples  on  applica¬ 
tion  of  the  Neutral  Action iMethod^ 


Acknowledgements 


c>r6z.r 


This  work  was  performed  with  the  support  of  the  U.S.  Department  of 
Energy  and  the  Air  Force  Office  of  Scientific  Research.  This  support  is 
gratefully  acknowledged. 


R<r«rtaca 

[1]  E.  Besiel-Hat^a.  Gbtr  die  ErhattunfuStze  der  Eltcirodyniunik.  Maih.  Ann.  S4,  158-276  (1921). 

(2)  E.  Noether,  liworianit  Vanciionsprobleme.  Nachr.  Kdnig.  Gea.  Wiiien.  GStiingen.  Matb.-Phys.  KJ. 
pp.  2}5-257(  1918)  (jeeTraruport  Theory  and  Sul.Phys.  I,  l86-207(]97l)ror  an  Engliih  translation). 

())  P.  J.  Giver,  Applications  of  Life  Croups  to  Differemust  Equations.  Springer,  New  York  1986. 
j«j  T.  Honeio,  N.  Chien  and  G.  Herrmann,  On  conservation  laws  for  dissipotiae  systems.  Phys.  Lett.  A. 
ISS,  223  -  224  (1991). 

(51  ht.  Ctuen,  T.  Honem  and  G.  Herrmann,  Dissipative  Systems,  Conservation  Laws  and  Symmetries,  1991, 
be  published. 

Ahttraet  , 

r,  O 

la  the  brief  noultitledl  “On  Conservation  Laws  for  Dissipative  Systems'*  141.  a  new  method  for 
eonstfueiine  conservation  laws  was  proposed.  This  method  wu  termed  tbef"NeutraJ  A  .Uon  tNA) 
Method  j  in  (5].  For  any  system  govern^  by  a  set  of  difTerenlial  equations,  the  NA  method  otters  a 
'lystematK  approach  for  deienninstion  o(  conservation  laws  applicable  to  the  system.  It  is  the  purpose 
of  the  present  paper  to  establish  conservauon  laws  for  one-  and  two-dimensional  viscoelasticy  (Voigt 
model)  via  the  NA  owthod.  The  conservation  laws  derived  should  prove  useful  in  studies  of  fracture  and 
defects  in  a  viacoelasiic  material. 
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